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Abstract 

The behaviour of the topological susceptibility % in QCD with two colours 
and 8 flavours of quarks is studied at nonzero temperature on the lattice across 
the finite density transition. It is shown that the signal of % drops abruptly at 
a critical chemical potential fj, c , much as it happens at the finite temperature 
and zero density transition. The Polyakov loop and the chiral condensate 
undergo their transitions at the same critical value \x c . At a value \i s of the 
chemical potential, called saturation point, which in our case satisfies fj, s > fj, c , 
Pauli blocking supervenes and consequently the theory becomes quenched. 
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I. INTRODUCTION 



The topological susceptibility x nas many applications in QCD. It is denned as the 
correlation function of two topological charge density operators Q(x) at zero momentum 
where Q(x) is 

Q{x) = ^c^W F W . (i-i) 

and F£ v (x) is the gluon tensor. The quenched value of x is responsible for the large mass of 
the r]' meson [1,2] and allows to understand the realization of the Ua(1) symmetry. From 
phenomenology this quantity must be x ~ (180 MeV) 4 . The most recent lattice results [3-5] 
and calculations based on SVZ-sum rules [6] are in fair agreement with this prediction. Be- 
sides, in the full theory x is proportional to the chiral condensate in the massless limit [7] 
while the second moment of the correlation function of two topological charge density oper- 
ators furnishes a possible explanation to the so-called proton spin crisis [8]. 

The behaviour of the topological susceptibility at finite temperature and/or matter den- 
sity is an important ingredient to understand whether or not the singlet chiral symmetry 
undergoes some sort of effective restoration at extreme conditions. The fate of this symmetry 
and of the non-singlet chiral symmetry and the ordering of the respective phase transitions 
open several scenarios [9-12] with relevant physically measurable consequences. 

It is known that x undergoes an abrupt drop at the deconfinement transition at finite 
temperature in Yang-Mills theory, both for two [13] and three [14] colours, and in full 
QCD for three colours and 2 or 4 flavours of staggered quarks [15]. These results have 
been confirmed by many collaborations by using cooling [3,16,17]* (see however [18]) and by 
counting fermionic zero modes [19], including studies for iV > 4 colours. 



* Older papers claimed that the signal dissapeared smoothly and it began to fade away well before 
the deconfinement temperature; however nowadays this possibility seems excluded. 
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In this paper we present an investigation about the behaviour of x & t finite temperature 
and density. This study has been carried out considering two colour QCD. In this case 
we do not have to face the hard sign problem of numerical simulations in three colour 
QCD at finite density (although several techniques to overcome these difficulties have been 
proposed [20-22]). 

Several analyses of two colour QCD at finite density are available in the literature [23-31] 
and one expects that the qualitative aspects of the behaviour of gluonic observables can be 
extended to three colour QCD. In particular, for the topological susceptibility, this expec- 
tation is further supported by the fact that instantons [32] are present in any SU(iV) gauge 
group theory. Monte Carlo results bear out this expectation both for the validity of the 
Witten-Veneziano mechanism as well as for the sudden drop of the signal of x & t the de- 
confinement transition [3,13,17]. 

Model studies suggest that instantons form structures with different features as the 
temperature and density change [33,34]. Beyond the finite density transition instantons 
seem to cluster into polymers while at the deconfinement temperature and zero density they 
tend to form instanton-antiinstanton pairs. 

We have investigated the behaviour of x on the lattice as the density of matter is varied at 
finite temperature. The temperature was introduced as usual by compactifying the temporal 
direction of the lattice and the density by switching on a chemical potential \i for quarks. 
Fermions have been discretized by using the staggered formulation [35] 

The two main results of our numerical study are: i) the topological susceptibility drops 
at a critical ji c and ii) this ji c is the same value of the chemical potential where the chiral 
condensate falls and the Polyakov loop rises. Moreover the introduction of a density of matter 
brings forth the following effect. The density of fermions grows as the value of [i is increased 
up to a saturation point, ji s > /i c , where no more fermions can be inserted into the system. 
In fact, due to the Pauli exclusion principle, there cannot be more than Nf/4 fermions per 
site for SU(2) gauge group and Nf flavours (which in our case is equivalent to Nf/8 baryons 
per site). This number follows from the fact that in the staggered formulation, one Dirac 
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spinor is reconstructed from the fields at the 16 sites of every hypercube [36] and that we are 
considering two colours. When the system is completely saturated with fermions, any further 
increase of the chemical potential will be totally ineffective: the fermions are completely 
frozen and the gluons are the only dynamical degrees of freedom. Then the observables 
take on their quenched values at /i > fj, s and, in particular, the topological susceptibility 
grows again until reaching a second plateau. In fact the corresponding temperature for the 
quenched theory lies in the confined phase where x is nonzero [14]. However it is important 
to notice that the saturation point is an unphysical effect related to the countable number 
of sites in the lattice or, stated otherwise, to the presence of an ultraviolet (UV) cutoff, and 
is thus expected to disappear in the continuum limit; an essential point is that in our case 
the saturation point /i s turns out to be larger than the critical chemical potential /i c and 
also well separated from it, so that the study of the deconfinement phase transition is not 
affected by that lattice artifact. 

The paper is organized as follows: in Section 2 the technical details of the numerical 
simulations are given; the calculation of x an d the investigation of its behaviour across 
the transition are presented in Section 3 together with the extraction of the physical units. 
Section 3 also contains a study of simultaneity of the chiral, confinement and topological 
transitions. Some conclusive comments are left for the last Section. 

II. MONTE CARLO SIMULATIONS 

We have performed numerical simulations in QCD with dynamical staggered quarks [35]. 
Following Ref. [23] we have worked with two colours and 8 flavours of staggered fermions. 
The action is [37] 

S = ^Y,Vi(x){e a ^(x)U 4 (x)^(x + Z) - e^ a ^(x + Z)Ui(xyip(x)\ 

1 3 _ 

+ 9 rii(x){il)(x)Ui(x)il)(x + i) - ip(x + i )C/ i (x)V(^)} 

^ i=l x 
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+am tp(x)tp(x) + Sw , (2.1) 

X 

where r] v (x) are the staggered phases and U v (x) are the gauge links starting at site x and 
directed towards V. Sw is the usual Wilson action for gauge fields [38]. Fermion fields 
have been multiplied by a suitable power of the lattice spacing a in order to render them 
dimensionless. 

For further details on the algorithm and the action we refer the reader to Ref. [23]. We 
remind the reader that the chiral symmetries of the lattice action are not equal to the ones 
of the corresponding continuum theory. 

We have numerically simulated the model on a 14 3 x 6 lattice at inverse bare gauge 
coupling P = 4/g 2 = 1.5 and bare quark mass am = 0.07 by using the Hybrid Molecular 
Dynamics (HMD) algorithm [39]. The presence of a chemical potential prevents us from 
separating odd and even sites [40] and as a consequence we worked with 8 flavours. We have 
not used the square root trick to reduce the number of quark species. 

We have run the Molecular Dynamics algorithm with step size St = 0.02 [23] and 50 
decorrelation HMD steps have been performed between two successive measurements. This 
separation was chosen so to have configurations with well-decorrelated topological proper- 
ties [41-43]. For a few values of //, 100 HMD steps were used. 

Although we have used a non-exact algorithm, possible systematic errors due to the 
finite step size seem to be under control. This statement is corroborated by two facts: i) 
runs corresponding to the same input parameters but with finer HMD step size (St = 0.005) 
have been carried out and results remain stable; ii) at large chemical potential, where the 
theory becomes effectively quenched, our results are in complete agreement with similar 
measurements obtained from separate runs performed by simulating the pure gauge action, 
(see Section 3.c). 

We have calculated the topological susceptibility from the fluctuations of the topological 
charge, estimated by using the 1-smeared operator (see next Section). We also measured the 
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FIG. 1. Xl? (circles) as a function of the adimensional chemical potential aji. The last 



(i) 



point (empty down triangle) is the value of x L 



at the same (3. In the inset the ratio 



quenched 

x(m)/x(a* = 0) (squares) is shown in the interval of a\x where no apparent effects of Pauli blocking 
are detectable. Lines are drawn to guide the eye. 



Polyakov loop, the average plaquette, the baryon density and the chiral condensate. De- 
pending on the value of a/i, the statistics for these observables was 300-500 uncorrelated 
measurements. 

Early exploratory runs indicate that on a 14 3 x 6 lattice at f3 — 1.5, am = 0.07 and zero 
chemical potential, the system is in the normal hadronic phase. According to the suggested 
phase diagram for two colour QCD in the temperature-density plane (Fig. 1 of Ref. [44]), 
our simulations will cross a transition line which separates the hadronic phase from another 
state. Depending on the temperature, this state can be either the quark-gluon plasma or 
the superfluid state and the respective transition can change its nature correspondingly. In 
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FIG. 2. Polyakov loop P as a function of afi. The logarithmic scale allows to disentangle the 
data obtained in the vicinity of the transition point. Points are joined by a line to guide the eye. 

the present work we have not examined the order of this transition but we have shown that 
the transitions for the Polyakov loop, chiral condensate and topological susceptibility are 
concomitant. From this fact and the calculated value of the ratio T/T c (see later) we infer 
that possibly our system becomes a deconfined plasma of quarks and gluons for \i > ji c (T c 
is defined as the temperature where the chiral vector symmetry is restored at zero chemical 
potential). Moreover the fermionic saturation is consistently observed at roughly the same 
value ji s for all observables. 
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III. MEASUREMENT OF X AND ITS BEHAVIOUR ACROSS THE TRANSITION 



A. The method 



The first task to calculate the topological susceptibility on the lattice is the introduction 
of a regularized topological charge density operator Ql(x). We require that this operator 
satisfies the continuum limit (in the operator sense) Q L{x)-—>a> A Q{x) where a is the lattice 
spacing. By calling Ql = J2xQl(x) the lattice topological susceptibility is defined as 



where V is the spacetime volume of the lattice. This quantity is not yet equal to the physical 
susceptibility X - They are related through the general expression [45] 



where Z and M are renormalization constants usually called multiplicative and additive 
respectively. They depend on the lagrangian (specifically on the quark masses and the 
gauge coupling) and on the operator used on the lattice to regularize the topological charge. 
In particular had we used an operator Ql(x) suggested by a fermion theory obeying the 
Ginsparg-Wilson condition [46,47] then Z would have become equal to 1 and M equal to 
zero [48] . We did not make use of such operators because their calculation on the lattice is 
very demanding in computer time. 

In order to obtain X from Eq. (3.2) we have used the Pisa method (also called "field- 
theoretical method"). In this method one first calculates Z and M and then inserts their 
values in Eq. (3.2), to extract the result for a 4 x- 

Let us explain the origin of Z and M in Eq. (3.2). Let G^ n ' l \pi,p2, • • • ; qi, qi • ■ •) be a 
bare Green's function with n elementary fields (carrying momenta pi) and I insertions of the 
topological charge density operator (carrying momenta qi) to be computed in the quenched 
theory. This function can be calculated either in the continuum by using some regularization 




(3.1) 



XL = Z 2 a 4 x + M, 



(3.2) 
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method or on the lattice. In both cases the inclusion of renormalization constants for the 
external fields and gauge coupling is enough to make the Green's function finite (because 
the calculation is performed in the quenched theory where Q(x) is a renormalization group 
invariant operator). Let us call G^^ttice an d C^continuum the two renormalized functions thus 
obtained. In general, for an arbitrary renormalization scheme, the two functions G^fl ttice 
and C^continuum are n °t necessarily equal and they match only after the inclusion of a mul- 
tiplicative finite renormalization constant Z [49,50] for each of the / insertions of topological 
charge operator. In a formal writing: Ql(x) = Za A Q(x). When I > 1 there may be further 
divergences originated by contact terms, see later. 

In the theory with fermions the topological charge mixes with other operators related 
to the axial anomaly [51]. This mixing induces a correction to the above described multi- 
plicative renormalization. Such a correction is however rather small [52] and we neglect it, 
consistently with the large statistical and systematic errors of our numerical simulation. 

The r.h.s. of Eq. (3.1) contains the product of two operators at the same spacetime 
point. In Field Theory such products can be divergent [53]. On the other hand the correlator 
(Q(x)Q(0)) is negative for x ^ as shown in [54-56]. Also on the lattice this correlator is 
negative [57-59]. Then part of the contact divergences must add up to the total topological 
susceptibility in order to make it positive. The rest of the contact divergences, if any, must 
be subtracted. This subtraction is M. In order to calculate it, we can follow the strategy 
introduced in [61,14]: M is the value of xl in the sector of zero total topological charge, 
M = Xl q q - This prescription guarantees the physical requirement that x vanishes in 
such a sector. M can be also calculated in nonzero topological charge sectors [62] yielding 
concordant results. Although M consists of UV divergences, it turns out to be a finite 
quantity because we are working on the regularized theory (in the practical calculation the 
cutoff a is never sent to zero). 

We have calculated the topological susceptibility for several values of the chemical po- 
tential at a fixed inverse gauge coupling (3. This entails that the additive and multiplicative 
renormalization constants must be the same for all values of //. Since we will present the res- 
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Table 1 



Measurements of several operators for varying \i. 



an 


io 5 xi 1} 


10 5 M 


a 3 p B 


(|TrD) 


(V^)full 


quenched 


P 


0.0 


7.77(95) 


2.54(19) 


0.001(1) 


0.4621(2) 


0.592(2) 


0.904(3) 


0.023(1) 


0.09 


7.57(49) 


— 


— 


0.4627(3) 


0.587(2) 


— 


0.0283(7) 


0.15 


7.46(48) 


— 


— 


0.4631(2) 


0.581(2) 


— 


0.0348(6) 


0.17 


7.60(49) 


— 


0.005(1) 


0.4647(4) 


0.568(2) 


— 


0.0472(9) 


0.18 


6.37(42) 


— 


0.006(2) 


0.4665(2) 


0.552(2) 


— 


0.0551(7) 


0.2 


4.24(47) 


2.29(19) 


— 


— 


— 


— 


— 


0.3 


2.88(21) 


— 


— 


— 


— 


— 


— 


0.4 


2.45(21) 


2.16(13) 


0.060(2) 


0.4912(2) 


0.280(1) 


0.493(3) 


0.1894(6) 


0.5 


3.54(25) 


— 


— 


— 


— 


— 


— 


0.6 


4.69(47) 


2.73(19) 


0.230(2) 


0.4764(2) 


0.105(1) 


— 


0.3085(6) 


0.7 


6.40(47) 


— 


— 


— 


— 


— 


— 


0.8 


8.50(60) 


— 


0.554(2) 


0.415(2) 


0.062(1) 


0.059(9) 


0.3289(7) 


1.0 


9.29(57) 


— 


0.882(2) 


0.3742(5) 


0.015(3) 


0.005(1) 


0.157(3) 


1.1 


— 


— 


0.958(1) 


0.3668(3) 


0.003(1) 


0.003(1) 


0.099(2) 


1.2 


9.59(62) 


— 


0.987(1) 


0.3638(1) 


0.0014(3) 


3(5) 10~ 4 


0.057(1) 


1.4 


9.98(67) 


— 


0.998(1) 


0.3626(4) 


0.0012(4) 


3(3) 10~ 4 


0.019(1) 


1.5 


— 


— 


0.999(1) 


0.3626(1) 


1(2) 10~ 4 


— 


0.010(1) 


1.6 


9.70(58) 


— 


1.000(1) 


0.3625(1) 


-2(2) 10~ 4 


-3(2) 10~ 4 


0.006(1) 


1.8 


9.03(59) 














2.0 


9.83(67) 




1.000(1) 


0.3625(1) 


3(5) 10~ 5 




7(2) 10~ 4 


2.4 






0.999(1) 


0.3624(1) 


6(6) 10~ 5 




1(4) 10" 4 


2.8 






1.000(1) 


0.3623(1) 


3(4) 10~ 5 




1(4) 10" 4 


3.2 






1.000(1) 


0.3626(1) 


2(3) 10~ 5 


-1(1) io- 4 


3(4) IO" 4 
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FIG. 3. Baryonic density a 3 ps (circles) and plaquette (^TrD) (up triangles) as a function of 
a/j,. The last point (empty down triangle) is the value of the plaquette for the quenched theory at 
the same (3 value. Lines are drawn to guide the eye. 

ults for the ratio x(A t )/x(A t = 0) as a function of //, the multiplicative Z needed not be 
calculated (the notation xil 1 ) indicates the physical topological susceptibility at a value \i 
of the chemical potential). Instead M was calculated by use of the heating method [60,61]. 
This method consists in the following: we started with the trivial configuration (all links 
equal to unity) which clearly belongs to the zero topological sector. Then we applied 5000 
steps of HMD updating and measured xl every 100 steps. This set of 50 measurements is 
called "trajectory". After each measurement we applied also 8 cooling steps [45,63,64] and 
determined Ql and the total energy of the configuration to verify that the topological sector 
was not changed (other methods to check the background topological sector, like the counting 
of fermionic zero modes [65], yield the same results for M and Z). We repeated the whole 
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. We estimated the 

Q=o 



procedure to obtain a number of trajectories. For each trajectory we always discarded the 
first few measurements because the configuration was not yet thermalized'''. Averaging over 
the thermalized steps (as long as the corresponding cooled configuration showed the correct 
topological charge sector, Q = within a deviation 5) yielded xl 
systematic error that stemmed from the choice of 5 as in Ref. [14] . 

When the cooling (or other) test reveals that a configuration has moved from the zero 
topological charge sector to another sector, it is discarded. Actually this event seldom 
happens because the topological modes are effectively decoupled from the UV modes so that 
after starting from a classical configuration of any fixed topological content, it is difficult to 
alter the background topological sector by applying some updating (heating) steps at the 
corresponding value of f3 to thermalize the UV fluctuations. In fact this decoupling of the 
two types of modes is the gist of the heating method [60,61]. 

B. The lattice operator 

A valid lattice regularization of the topological charge density operator, Eq. (1.1), is 
given by [66] 

1 ±4 

Q L W = -$r^ H e (uia Tr {U (u {x)U ya {x)} , (3.3) 

<^7<r=±l 

where H^ v (x) is the plaquette in the ( — v plane with the four corners at x, x + (, x + ( + V, 
x + V (counter-clockwise path). As indicated in the sum, indices (, u, 7 and a can point 
towards either positive or negative directions. Links pointing to negative directions mean 
U-„(x) = Ul(x — V). The generalized completely antisymmetric tensor is defined by €1234 = 1 

and C^—Qv^a = t^vya- 



tin this context thermalization is defined with respect to short range fluctuations within the zero 
topological charge sector. 
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1 




FIG. 4. Chiral condensate for the full theory (circles) and the quenched theory (triangles) as 
a function of afi. The logarithmic scale allows to disentangle the data obtained in the vicinity of 
the transition point. Lines are drawn to guide the eye. 

The lattice operator used in our simulations was the 1-smeared Ql\x) which is con- 
structed from Eq. (3.3) after substituting all link matrices by 1-smeared links [67]. If U v {x) is 
the link starting at site x and pointing towards the v direction, the 1-smeared link f/W [x) 
is denned by the following linear combination of the seven shortest paths connecting the 
points x and x + 9 

c ±4 

U v (x) = {l-c)U v {x) + - ]T U a (x)U„(x + a)Ul(x + V), 

6 Q t±l 
\a\& 

P 'W (3.4) 
(lTrUUx)U„(x)) 

This combination is then projected back onto SU(2) as indicated by the second line in 

13 



Eq. (3.4). By using smeared operators less noisy results are obtained. The parameter c can 
be tuned in order to optimize this improvement. We chose c = 0.85 as in [13]. The naive 
continuum limit of Ql\x) equals the continuum topological charge density, Eq. (1.1). 

Following the definition (3.1) we called = J2 x Ql\ x ) the total 1-smeared lattice 
topological charge and 

the 1-smeared lattice topological susceptibility. Xl? an d fl4 X are related by Eq. (3.2). We 
have calculated M for four values of \x as shown in Table 1. Averaging over all four results 
leads to M((3 = 1.5, am = 0.07) = 2.37(8) x 10~ 5 . This is the value of M that we have used 
throughout the paper. 



C. The Monte Carlo simulation 

In Fig. 1 the results for x*l are shown as a function of a/i. Data are given in Table 1. 
The sudden drop at a/j, c = 0.175(5) is apparent. In the inset of Fig. 1 we show the scaling 
ratio x(A t )/x(/ u — 0) f° r < 0-4- Accordingly with Eq. (3.2) this ratio is calculated as 

X(A* = 0) xi 1} (^ = 0)-M' 
The notation Xl\^) stands for the lattice topological susceptibility, Eq. (3.5), calculated at 
chemical potential \i. This inset is one of the two main results of our paper: it displays the 
drop of the physical topological susceptibility. Notice that the results in this plot are free 
from Z and a 4 factors. 

In Fig. 2 the Polyakov loop is shown in a similar range of values of chemical potentials 
at the same f3 and quark mass. All data are listed in Table 1. The rise of the Polyakov loop 
displayed by this Figure happens at a value of a/j, which is compatible within errors with 
the value of a\i c at which xip) drops. 

For values of the chemical potential larger than a/i = 0.5, the topological susceptibility 
begins to increase again and then flattens at afj, > 1.2, as shown in Fig. 1. Such a value 
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corresponds to the saturation point, ji s . This new plateau is not the indication of a new 
physical phase. It is due to the Pauli blocking, a saturation of the fermion number through- 
out the lattice. The value for the chemical potential is so high that the baryon density 
Pb reaches the maximum allowed value, 1 baryon per site. After this point the fermions 
are frozen and the theory is entirely analogous to the quenched theory. To say it plainly, 
simulating the theory at those high values of /i is nothing but a clumsy way to reproduce 
the SU(2) pure gauge theory. 

The above qualitative argument is validated by Monte Carlo simulations of the pure 
gauge theory that compare very well with the full theory simulated at /i ^ /i s . In particular 
the empty down triangle in Fig. 1 is the value of x^l calculated in a separate simulation by 
using the SU(2) pure gauge Wilson action, its value being x [ l\ = 9.84(19) 10" 5 . The 

quenched 

agreement between this number and the last data for a/x ^ 1.2 is manifest. Also the signal 
for the Polyakov loop in Fig. 2 disappears at a/i > 1.2. This is the expected behaviour 
for this operator in the quenched theory when, as in our case, the temperature is below 
deconfmement, T < T c . 

The rise of the signal of Xl^ f° r I 2 ~ i s s l° w i n contrast with its sudden drop at /i c . 
This effect is due to the partial blocking of fermions which increases with \i until it reaches 
the total saturation at ji s . 

Notice that the saturation of the fermion degrees of freedom makes the theory to behave 
as the pure gauge theory in all respects. In particular not only Xl? assumes the quenched 
value, but also the renormalization constants and the beta function do (i.e.: the full theory 
becomes quenched, also perturbatively), in such a way that the physical susceptibility x 
turns into that of the pure gauge theory too. 

The saturation of fermions and subsequent forced quenching of the theory is due to the 
fact that on the lattice there are a finite number of physically accessible sites. Hence it 
is an unphysical lattice effect with no counterpart in the continuum. It is interesting to 
notice that the effective quenching at large values of // > /j, s corresponds to the fermionic 
determinant becoming a constant independent of the gauge field background; this can be 
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easily understood looking at the fermionic part in Eq. (2.1): in the large p limit the forward 
temporal part of the Dirac matrix is exponentially larger than the other terms, and the 
determinant of a matrix containing only the forward temporal part is indeed independent 
of the gauge field background, as can be easily verified. 

In Fig. 3 we show the data for the baryonic density a 3 p B and the plaquette (|TrD) versus 
the chemical potential. Both data sets can be found in Table 1. At the saturation point, 
a/i s ~ 1.2 , the density attains the maximum value a 3 pB = 1 and the plaquette becomes 
that of the quenched theory. This value of a/i s is roughly the same at which the topological 
susceptibility begins the second plateau in Fig. 1 and the Polyakov loop becomes zero in 
Fig. 2. The quenched value of the plaquette at f3 — 1.5 has been calculated in a separate 
simulation with the pure gauge Wilson action and it is shown as an empty triangle in Fig. 3. 
It amounts to (|TrD) quenc h e d = 0.362432(23) which is in perfect agreement with the (Pauli 
blocked) values of the average plaquette in the full theory for /x > p s , (see Table 1). 

All observables agree with their corresponding quenched values as soon as the chemical 
potential is high enough to induce fermionic saturation. We have also calculated the chiral 
condensate as a function of ap. The results are shown in Table 1 and Fig. 4 for both the full 
and quenched theories. In the two simulations, we have measured the same operator, i.e. the 
trace of the inverse of the (/i-dependent) quark matrix. In one case it was calculated by using 
configurations obtained from the full theory and in the other case by using configurations 
corresponding to the pure gauge theory. In the full theory the condensate drops at the same 
ap c where the Polyakov loop rises and the topological susceptibility falls abruptly. Thus the 
three operators undergo their transitions at the same density. We are aware that neither the 
Polyakov loop nor the topological susceptibility are order parameters for confinement and 
the singlet flavour Ua(1) symmetry respectively. In particular the drop of the signal of the 
topological susceptibility does not necessarily mean a restoration of the Ua(1) symmetry 
although it must have physically measurable consequences in observables directly related 
to x- This set of consequences is commonly termed as "effective restoration" of the Ua(1) 
symmetry [9,11]. 
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For values of the chemical potential larger than the saturation point /x s the chiral con- 
densate in the full theory stays vanishingly small. This behaviour coincides with that of the 
chiral condensate calculated in the quenched theory as shown in Fig. 4. Notice that in the 
present case, the operator depends on the chemical potential and this is why we show the 
data for the quenched model as a function of a/i. However as soon as the chemical potential 
reaches the saturation value the two sets of data become coincident. This is another indi- 
cation that the full theory at /i > /i s becomes entirely quenched. In the quenched case a 
straightforward algebra yields the leading dependence for large \x 

(^quenched OC e~ Nt ^ , (3.7) 

where N t is the lattice temporal extent, N t = 6 in our case. 



D. The numerical results 

We have also given approximate physical units to the lattice spacing a and therefore also 
to the critical value of the chemical potential /i c . In a separate Monte Carlo simulation we 
have calculated the chiral condensate and the Polyakov loop as a function of the temperature 
at zero chemical potential. To this end we have kept fixed the temporal lattice size and have 
varied the gauge coupling (3. In Fig. 5 and Table 2 we show the results of the simulation. 
From the position of the flex point in any of the two curves of this Figure we can extract 
the critical beta, (3 C = 1.594(6). By taking into account the two loop beta function of the 
model (Al is the renormalization group invariant mass scale corresponding to the bare lattice 
coupling of the action (2.1)), 

a (^ = e "* a "fe) 5/2 [ 1+0 G)]' < 3 - 8 > 

we obtain the ratio a((3 c )/a(/3 = 1.5) = 0.34(3) where the error comes from the imprecision 
in the determination of j3 c (notice that we use indistinctly the notation a or a((3)). The error 
caused by the ignorance of higher loop terms in the beta function is rather small within the 
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narrow interval between (3 — 1.5 and (3 C (actually the sole one loop term is enough to account 
for about the 85% of the variability of the function a((3) in that interval). 

Since our simulations have been performed at (5 — 1.5 on a lattice of time size equal to 6, 
the temperature of our system was T/T c = a((3 c )/a((3 = 1.5) = 0.34(3) and the resulting 
value for the critical chemical potential was p c /T c = 6 x (T/T c ) x 0.175(5) = 0.357(10)(32) 
where the first error is due to the Monte Carlo determination of the critical value of ap and 
the second one derives from the error in j3 c . 

From Table 1 we can also give an approximation for the critical density. Since the 
transition lies between ap = 0.17 and ap = 0.18, from Table 1 we can approximate the 
density at the transition point as a 3 p Bc = 0.0055(15) which is the average of a 3 ps among 
the two values of ap. From this average we obtain Pb,c/T 3 = 6 3 x (T/T c ) 3 x 0.0055(15) = 
0.047(13) (12) where the first error comes from the lattice determination of a 3 ps c and the 
second one is due to the determination of f3 c . 

Based on the large number of fermions in our lagrangian and on the fact that they 
are quite heavy, we surmise that the critical temperature for our system at zero chemical 
potential is T c = 100—200 MeV. Then the value of the lattice spacing at f3 c is a((3 c ) = 1/6T C = 
0.00111(^8) MeV" 1 . From the ratio T/T c we obtain the value of the lattice spacing m our 
simulation runs, a((3 = 1.5) = 0.00327(22) (t^ 5 ) MeV" 1 = 0.64(4) (t?jj) fm where the first 
error is due to the imprecision on (3 C and the second one to the estimated inaccuracy on T c . 

We see that our lattice is coarse. It is a consequence of the small lattice size (due to 
computer limitations) and of the need to work at T < T c in order to run into a phase 
transition while p is increasing. 

Notice that within the above limits for T c , the value of the critical fermion chemical 
potential p c ranges from 36(1) (3) to 72(2) (6) MeV. Then the critical baryonic chemical 
potential p B , c — 2/i c varies from 70 to 145 MeV at a temperature T/T c = 0.34(3), indicating 
a curvature of the phase boundary in the p-T plane quite larger than what is measured 
in QCD with three colours: that does not come as a surprise, since in two colour QCD 
baryon states are degenerate with mesons [68] , so that the critical chemical potential at zero 
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temperature may well be lower than in usual QCD. 



Table 2 

Measurements of (^^)m\ and P for varying (3 at ji — 0. 





(V^)full 


P 


1.525 


0.5521(19) 


0.030(1) 


1.55 


0.5080(25) 


0.038(1) 


1.5625 


0.4828(25) 


0.0402(35) 


1.575 


0.4466(26) 


0.0549(23) 


1.5875 


0.4131(30) 


0.0685(10) 


1.6 


0.3735(22) 


0.0880(18) 


1.6125 


0.3454(24) 


0.1041(20) 


1.625 


0.3279(15) 


0.1082(20) 


1.65 


0.3005(8) 


0.1208(10) 


1.675 


0.2801(7) 


0.1315(7) 



By allowing the valuation of T c to vary in such a wide range (100-200 MeV) we have 
mimicked an estimate of the systematic error derived from the great number of quark species 
and their big mass. There are however other sources of systematic errors which have not 
been contemplated in the present analysis: the (wrong) gauge group, the large value of the 
lattice spacing, the modest size of the lattice volume and the use of a non-exact updating 
algorithm. The experience in simulating and comparing the physics of SU(2) and SU(3) 
gauge theories teaches us that the first problem is possibly unimportant with regard to the 
behaviour of the topological susceptibility. As described above, the main influence of the 
wrong gauge group regards the numerical values of the physical parameters. The second and 
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FIG. 5. Chiral condensate (circles and left vertical axis) and Polyakov loop (squares and right 
vertical axis) as a function of (3 at vanishing chemical potential in the neighbourhood of the finite 
temperature transition. The vertical line and band are the value of (5 C and its error. 

third problems can be ameliorated by studying the theory on larger volumes. We think that 
none of these errors can change the two main conclusions of our paper: the drop of the 
topological susceptibility at the transition and the coincidence of all transitions. 



E. Relationship between transitions 

We have claimed that the three transitions studied in the present paper, namely the 
one related to topology, to the Polyakov loop and to the chiral condensate are concomitant: 
they all happen at about the same value of the chemical potential. We substantiate this 
statement with the help of the plot shown in Fig. 6. In this plot we display the derivatives 
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0.1 0.2 0.3 0.4 

FIG. 6. Derivatives with respect to the chemical potential jjl of the normalized topological 
susceptibility (circles, continuous line and left vertical axis), Polyakov loop (up triangles, dotted 
line and right vertical axis) and chiral condensate (down triangles, dashed line and right vertical 
axis). The two vertical axes are expressed in units of the lattice spacing a. The three sets of data 
have been slightly shifted horizontally to avoid the overlapping of various symbols and error bars. 
Lines are the result of the interpolation described in the text. 

of the three observables (with immaterial proportionality constants) with respect to the 
chemical potential. If the transitions are identified as the points where the variations are 
sharper then the Figure seems to give unquestionably evidence in favour of the coincidence 
for the three transitions. This conclusion is the second main result of the paper. 

The plot is shown within the interval of chemical potentials spanning from to the 
maximum \i where no Pauli blocking effects are apparent. The isolated points (circles, up 
and down triangles) for the observable O (x(/f)/x(0), P or (V'V')fuii) have been obtained by 
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calculating the discrete derivative 

CHjM) - 0(^_i) ^ g ^ 

— l^i-l 

from the data in Table 1 and placing the result at position a(/ij + /ij_i)/2. Furthermore 
Monte Carlo data for each observable have been interpolated by a natural cubic spline [69]. 
Then the derivative of these interpolation functions have been calculated and they are the 
lines plotted in Fig. 6. All three lines and all three sets of data points show a clear peak 
at the same value a/i c = 0.175. The continuous line, which corresponds to the topological 
susceptibility, displays other points with a large derivative. This is an spurious effect on the 
spline interpolation due to the fluctuation that makes the value of xl at a/i = 0.17 to move 
upward a little bit, see Fig. 1. 

A similar figure can be drawn for pure gauge theories and for full QCD indicating that 
the concomitancy is a common feature of the three kinds of transition for all theories that 
we have studied up to date [70]. 

IV. CONCLUSIONS AND OUTLOOK 

We have studied on the lattice the behaviour of the topological susceptibility x i n QCD 
with two colours and 8 flavours of quarks as a function of the quark chemical potential jj, at 
finite temperature. The so-called sign problem hinders a similar investigation in full QCD 
with three colours. 

Previous studies of topology in Yang-Mills theories at finite temperature for two and 
three colours both with or without fermions have shown that x undergoes an abrupt drop at 
the deconfmement phase transition. The results of our paper indicate that a similar quick 
drop takes place also at the finite density phase transition, (see inset in Fig. 1). It happens 
at a value of the chemical potential afx c = 0.175(5). This value of \i is the same at which the 
Polyakov loop and the chiral condensate undergo sudden changes (see Fig. 2 and Fig. 4). 
The rise of the Polyakov loop demonstrates enhanced screening possibly compatible with 
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deconfinement and the chiral condensate exhibits a behaviour that is compatible with a 
restoration of flavour chiral symmetry (see Fig. 4). 

If the chemical potential is further increased, the topological susceptibility rises again as 
more and more fermions are inserted into the system. We have given compelling arguments 
showing that this rise is not the signal of a new physical phase but just an artifact of the 
UV regularization of the lattice, which makes the number of accessible fermion levels finite. 
Indeed at the saturation point ji s (for our choice of parameters aji s ~ 1.2) the system is 
completely packed with fermions whose dynamics, due to the Pauli principle, gets frozen. 
Therefore for jj, > fj, s the gluons are the only dynamical degrees of freedom and the system 
behaves as the pure gauge theory at the same value of the gauge coupling and consequently 
the topological susceptibility takes the corresponding quenched value. In order to support 
these assertions, we have calculated the baryonic density and verified that it becomes equal 
to 1 baryon per site for all /i > /j, s , (see Fig. 3). Similarly, also other observables like the 
Polyakov loop (Fig. 2), the average plaquette (Fig. 3) and the chiral condensate (Fig. 4) 
take their quenched values for /i> /i s . 

We stress that we defined /x s as the value of the quark chemical potential where saturation 
becomes complete or, stated otherwise, where the theory becomes quenched. 

Since the Pauli saturation of fermions is a pure lattice effect with no counterpart in the 
continuum, we have not made a detailed study to establish whether its onset corresponds to a 
true (unphysical) phase transition or not: that may be the subject of further investigation. In 
fact n s must diverge in the continuum limit. Unfortunately this lattice artifact precludes us 
from making contact with the physics of extremely high densities by numerical simulations. 

We have shown that a drop in the signal of the topological susceptibility occurs in the 
theory with two colours and 8 equal mass flavours. What about the true theory? Due to 
the similarities in the topology of SU(2) and SU(3) it is conceivable to expect that the drop 
must be a genuine fact of QCD. Instantons are responsible for the topological properties 
and they are present in all SU(iV) gauge models. 

We call T c the temperature where the chiral vector symmetry is restored at vanishing 
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chemical potential. By resorting to the above-mentioned similarities between the present 
model and the true theory of the strong interactions, we have assumed that T c ranges 
between 100 and 200 MeV. From this approximation and knowing the value f3 c of the gauge 
coupling at which the chiral symmetry is restored for \i = (see Fig. 5), we have extracted 
an estimate of the lattice spacing in physical units. Also the ratio /j, c /T c = 0.357(10) (32) has 
been calculated, where the first error is due to the lattice determination of a\x at the critical 
point and the second one is derived from the error band in f3 c (see Fig. 5). Furthermore 
T/T c = 0.34(3). Notice that this ratio is insensitive to the imprecision on T c . Moreover 
the baryon density ps at the critical point is ps jC /T c 3 = 0.047(13) (12) where the errors 
have analogous meanings than those for fi c /T c : lattice evaluation of p^.c and error on f3 c 
respectively. All the above numerical results are subject to several sources of systematic 
errors (wrong gauge group, small lattice size, large lattice spacing). Hence we consider them 
as rough estimates of the true values. However none of the preceding sources of systematic 
errors is expected to be able to influence the main conclusions of the present paper, namely 
the drop in the topological susceptibility and the coincidence of this drop with the chiral 
and Polyakov loop transitions. 

By using the assumed limits to the value of the critical temperature at zero chemical 
potential T c , the critical baryonic chemical potential /ib,c = 2/z c can be estimated to range 
from 70 to 145 MeV at T/T c = 0.34(3). Therefore the value of pi B , c at zero temperature 
may be quite lower than what is expected in QCD with three colours. That is not a surprise 
since in the theory with two colours the lightest baryon state is degenerate with the pion. 

We have not studied the order of the transition and the properties of the states at both 
sides of this transition. On the other hand the number and mass of the quarks in our study 
are unphysically large. These are important aspects which deserve further analysis because 
the order of the transition might depend on these parameters. 

Our results indicate that the screening of colour forces, chiral symmetry and topology 
are interrelated in the T m (0.3 — 0.4)T C region of the two colour QCD phase diagram [44]. 
It would be interesting to repeat this study at other temperatures, in particular at very 
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small ones where the transition goes from the hadronic to the superfluid phase, which is 
still confining [71,72], and where chiral symmetry remains broken by a diquark condensate. 
Such low temperatures, say T < 0.1 T c , were not accessible to our simulations because on 
our lattice sizes it would have required such a large value of the lattice spacing that the 
results would probably have lost all physical meaning. 
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